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ABSTRACT 

High-resolution (R = 143,000), high signal-to-noise (S/N = 700-1100) 
Gemini-S bHROS spectra have been analyzed in a search for 6 Li in 5 stars which 
host extrasolar planets. The presence of detectable amounts of 6 Li in these ma- 
ture, solar-type stars is a good monitor of accretion of planetary disk material, 
or solid bodies themselves, into the outer layers of the parent stars. Detailed 
profile-fitting of the Li I resonance doublet at A6707.8 A reveals no detectable 
amounts of 6 Li in any star in our sample. The list of stars analyzed includes 
HD 82943 for which 6 Li has been previouly detected at the level of 6 Li/ 7 Li = 
0.05 ± 0.02. The typical limits in the derived isotopic fraction are 6 Li/ 7 Li < 
0.00-0.02. These upper limits constrain the amount of accreted material to less 
than ~ 0.02 to 0.5 Jovian masses. The presence of detectable amounts of 6 Li 
would manifest itself as a red asymmetry in the Li I line-profile and the derived 
upper limits on such asymmetries are discussed in light of three-dimensional hy- 
drodynamic model atmospheres, where convective motions also give rise to slight 
red asymmetries in line profiles. 
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( 119971 ): stars with giant planets are systematically metal-rich (by ~ 0.2 dex) relative to field 
FGK dwarfs not known to harbor planets. Two hypotheses have been proposed to explain 
this excess: primordial enrichment or pollution. The first indicates that the probability 
of forming giant planets is a steeply rising function of the intrinsic metallicity of the gas 
and dust cloud which gave birth to the system . This possibility is in agreement with the 
core-accretion scenario (e.g.. iPollack et al.lll996l ). The higher metal content would raise the 
surface density of solid material in the disk, leading to a more efficient agglutination of the 
cores onto which the gas will be accreted. The pollution scenario, on the other hand, indi- 
cates that the presence of planets could alter the metallicity of their hosting stars. During 
the inward migration process of giant planets, solid material from the protoplanetary disk 
or even inner planetesimals and planets could be accreted into the convective envelope of 
the hosting star. As this material is depleted in H and He, the star's metallicity would be 
enhanced. 



Current results (e.g., iFischer fc Valentil 120051 ) show that the frequency of planets in- 
creases significantly for higher metallicities, thus giving strong support for the primor- 
dial hypothesis. Evid ence for the occurence of pollution is still ambiguous. For instance, 
Ecuvillon et al.l (120061 ) studied the relation between chemical abundances of several elements 
and their respective condensation temperatures in stars with and without pl a nets, finding no 
signifcant differences in the two groups. On the other hand, iPasquini et al.l (120071 ) analyzed 
the metallicity distributions of planet-hosting dwarfs and giants and found that the latter do 
not favor metal-rich systems. The authors argue that this result could be a strong indication 
of pollution, as the metal excess could be erased by the dilution process that takes place 
during the later stages of stellar evolution. 

Another ambiguous result is the possible detection of 6 Li in the atmospheres of stars 
with planets, which is a sensitive test of the pollution hypothesis. Both lithium isotopes are 
destroyed at relatively low temperatures (T = 2 x 10 6 K for 6 Li and T = 2.5 x 10 6 K for 7 Li) 
in stellar interiors. During the early stages of evolution in solar-type stars (specifically, before 
entering the main-sequence), these stars are entirely convective and most of the primordial 
Li is transported to deeper and hotter layers, where it is rapidly burned. The fraction of 
lithium destruction is, however, a strong function of the stellar mass. For a given metallicity, 
there is a mass range in which the 6 Li is completely destroyed, while a significant amount of 
7 Li is preserved. The lower edge of this range falls just above the solar mass (corresponding 
to main-sequence late-F spectral types) for near solar metallicities. Thus, one should not 
expect to find any 6 Li in the atmospheres of solar-type stars. Any positive detection could 
indicate an external contamination or pollution process. 



Israelian et al.l (120011 ) measured the isotopic ratio 6 Li/ 7 Li in HD 82943 (with two close-in 
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giant planets) and found 6 Li/ 7 Li = 0.126 ± 0.014. This can be compared to a solar system 
(meteoritic) ratio of 6 Li/ 7 Li = 0.08. The positive 6 Li detection for th is star was interpreted 
as observational evidence of the pollution process. iReddy et al.l (120021 ) studied 6 Li in 8 planet 
hosting stars (HD 82943 included) and found no significant amount of this isotope in HD 
82943, nor in any of the targets analyzed. The difference from the previously published 6 Li 
detection for HD 82943 was attributed to the use of a more compl ete line list, al t hough 
it noted the presence of an unidentified absorption in the Li region. Ilsraelian et al.l (120031 ) 
investigated the nature of the unindentified absorption feature at 6708.025 A, which affects 
the Li I line, by observing several stars of differ ent effective t e mper atures. They concluded 
that a high excitation Si I line first propo sed by iMiiller et al.l (119751 ) is more adequate than 
the Ti I line used by IReddy et al.l (120021 ). Adopting a revised line list and higher quality 
spectra, the authors performed a new analysis of HD 82943 and measured 6 Li/ 7 Li = 0.05 ± 
0.02. Thus, this most recent result for HD 82943 gave additional support to the pollution 



scenario. 



More recently, iMandell et al.l (120041 ) made a major extension of the previous line lists 
(especially for the CN contribution) and tested three different possibilities for the unidentified 
feature at 6708.025 A (Si I, Ti I and Ti II): for all of the 3 different line lists no 6 Li was 
detected in a sampl e of three planet-hosting stars. These results are generally consistent with 
Reddy et al.l (120021 ) and argue against the pollution scenario. Unfortunatelly, HD 82943 was 
not analyzed by Mandell et al. (2004). 

In this paper, the 6 Li/ 7 Li isotopic ratio in HD 82943 is analyzed, as well as in other four 
planet hosting stars and one star not known to have giant planets. The observations and 
the data reduction are described in SJ2J The analysis procedures including the derivation of 
atmospheric and broadening parameters and compilation of the line list for the Li I region 
are presented in §31 The abundance results are presented in £@] and discussed in §5] 



2. OBSERVATIONS AND DATA REDUCTION 
2.1. Observations 

Spectra of the program stars were obtained at the Gemini-S telescope with the bench- 
mounted High- Resolution Optical Spectrograph (bHROS). Given the brightness of these 
targets, the "object-only" mode was used for observations; for this mode, a 0.9" fiber is 
fed into an image sheer that produces a 'slit', measuring 0.14" X 6.5" as projected to the 
camera focal plane. The spectrograph is cross dispersed by a set of fused silica prisms and 
an image sheer rotation mechanism is used to produce a 'vertical' slit at the observed central 
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wavelength on the detector, a single 2048 X 4608 E2V CCD with 13.5 /im pixels. 

The instrument was configured to produce a central wavelength on the detector of 6501 
A. In this configuration, eight incomplete spectral orders were obtained, covering the interval 
between 5580-7230 A. The continuous coverage available within a single echelle order varies 
from ~ 50-70 A depending on the order. The detector was used with lxl binning to achieve 
a resolution close to the spectrograph's nominal resolution of R = X/5X = 150,000 (3 pixel 
sampling) and minimize the impact of cosmic rays. The actual resolution was measured using 
10 ThAr lines in the order containing the Li I feature and this resolution was found to be R 
= 143,000 ± 5000. The spectra were obtained in a wide range of observing conditions during 
both classical (2006 May) and queue operations (2006 December and 2007 January). Tabled] 
contains a detailed log of the observations including spectral types, V magnitudes, number of 
exposures, integration times, and the resulting signal-to-noise ratios (per resolution element). 
The quoted S/N values are based on direct measurements of the rms in sample continuum 
regions which were selected based on inspection of the Solar Atlas. The measured values of 
signal-to-noise are compatible with the S/N based on Poisson statistics. 

In addition to the target spectra, calibration sets of biases, flats, and ThAr arcs were 
obtained each night. When the observations were done in queue mode, only a single ThAr 
spectrum was taken each night as only a single object was observed per night. When observa- 
tions were done in classical mode, ThAr spectra were taken which bracketed the observations. 
For example, the observations of HD 82943 obtained on May 9 2006 were bracketed with 
ThAr spectra and this provides a measure of the spectrograph stability. A direct measure- 
ment of the stability can be found when combining multiple spectra via crosscorrelation. The 
largest drift found was about 3 mA or 1/5 of a pixel. Telluric spectra were also obtained 
and, as expected, demonstrated no significant atmospheric contamination in the region of 
the lithium feature. 



2.2. Reduction 

The data were processed using standard echelle reduction practices with IRAP^ pack- 
ages, following the recipe developed by the bHROS science demonstration teauB The raw 
spectra were corrected by overscan and bias-image subtraction. The images were then flat- 



1 IRAF (Image Reduction and Analysis Facility) is distributed by the National Optical Astronomy Ob- 
servatories (NO AO), which is operated by the Association of Universities for Research in Astronomy, Inc. 
(AURA) under cooperative agreement with the National Science Foundation (NSO). 

2 Available at |http://www.gemini.edu/sciops/instruments/hros/hrosIndex.html| 
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fielded to remove both pixel-to-pixel variations and small amounts of fringing present in the 
spectra; since the illumination pattern on the detector is nearly identical for the flat and tar- 
get spectra, the fringe pattern can be very well removed. Scattered light was then sampled, 
fit, and subtracted from our images before extraction. 

The use of the prism cross disperser combined with the long 'slit' length causes the 
spectral orders to become tilted away from the central order. The severity of this tilt can 
potentially result in the loss of resolution in the extraction process. To reduce the impact of 
the tilts, each of the eight spectral orders was divided into 14 subapertures, for a total of 112 
resulting subapertures. Each subaperture was individually extracted with optimal extraction 
method and blaze corrected with similarly extracted flat-field spectra. The ThAr spectra 
were also subdivided into the same subapertures and wavelength solutions were derived 
for each of the 112 subapertures, with a typical rms scatter of 0.001 A. These wavelength 
solutions were applied to the target subapertures and combined to produce a reduced one- 
dimensional echelle spectrum composed of eight spectral orders. Rejection methods employed 
in the combination process ensured the removal of cosmic rays from the final spectra. Finally, 
the spectra were continuum normalized and corrected for ra dial velocity s h ifts. M ore details 



about the bHROS data reduction process can be found in ISchuler et al.l (120081 ) . The final 



processed spectrum for a target star is shown in Figure [U as an example. 



3. ANALYSIS 

The first step in the analysis is to derive effective temperatures, surface gravities, and 
microturbulent velocities for the sample stars. In addition, broadening parameters affecting 
the observed spectral lines also need to be defined. A crucial point in the dertermination of 
6 Li/ 7 Li isotopic ratios is the construction of a detailed line list for the region around the Li I 
feature. This section presents a discussion of the analysis method adopted in this study. 



3.1. Line Selection and Derivation of Stellar Parameters 



The bHROS spectra have incomplete wavelength coverage (Section I2.ip which restricts 
the selection of iron lines to derive stellar pa rameters and me tallicities. A sample of Fe I 
and Fe II lines was compiled fro m the list in iTheveninl (|1990l ) and equivalent widths were 



measured in the solar spectrum (IKurucz et ajj 



19841 ) in order to select suitable lines and 



^/-values which produced an abundance scatter of less than or equal to 0.05 dex. The final 
list of adopted Fe lines, as well as the measured equivalent widths for the target stars and 
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Fig. 1.— The Gemini bHROS spectrum of HD 82943 (R ~ 143,000 and S/N ~ 1130), 
showing part of the spectral order which contains the lithium feature. The main spectral 
lines in this region are identified. 
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the Sun are presented in Table [2j The wavelenghts, lower excitation po tentials (LEP), an d 
(//-values were taken from the Vienna Atomic Line Database^ (VALD-2; Kupka et al. 1999 ). 
As a reference note the adopted line list yields a solar abundance A(Fe£F 7.48 ± 0.05 and 
£ = 1.24 km s -1 using a Kurucz model atmosphere with T c g- = 5777 K, logg = 4.44, £ = 
2.00 km b _1 , and \/E p = 1.25. 

Stellar parameters for the target stars were derived spectroscopically and followed stan- 
dard techniques. Effective temperatures were obtained from zero slopes in diagrams of Fe 
abundance versus excitation potential and surface gravities from ionization equilibrium be- 
tween Fe I and Fe II species. The microturbulent velocities were varied until the slopes of 
A(Fe I) versus log(WA/A) were zero. The abundan ces were derive d in LTE using an updated 
version of the spectrum synthesis code MOOg| ( ISnedenl Il973l ). The model atmospheres 
adopted in the analysis were interpolated from the ODFNEW grid of ATLAS9 models^ 
( iCastelli fe Kuruca 12004 ) . Figure [2] shows the iterated result of A(Fe I) versus LEP (top 
panel) and A(Fe I) versus \og(W\/\) (bottom panel) for the star HD 82943, as an example. 
The adopted stellar parameters for all target stars are given in Table [31 Our adopted values 
for effective temperatures are in good agreeme nt (within ^100 K) w i th the derived values 



using the photometric correlations presented in lRarmrez fe Melendezl (120041 ) . Differences of 
this amount have no significant effect in the derived isotopic fractions. 

The Fe abundance uncertainties can be estimated from the errors in T e g, logg, £, line 
gf- values, continuum placement, and equivalent width measurements, among others. We 
estimate the uncertainties in T c g to be ±100 K; log g ± 0.2 dex; and £ ± 0.2 km/s. The 
combined uncertainties in the derived Fe abundances are typical for this type of classical 
spectroscopic analysis and are estimated to be about ±0.1 dex. 

The luminosities of the program stars were calculated using apparent V-magnitudes, 
Hipparcos parallaxes, bolometric corrections from Giradi et al. (2002), and using M bol 
= 4.77. No interstellar extinction was included as the most distant star has d = 65 pc. 
Utilizing these luminosities and the T e g- values from Table [3J stellar masses were derived by 
placing the stars in a grid of evolutionary tracks from Girardi et al. (2000) for [Fe/H] = 
0.00 and 0.20, with stellar masses covering M = 0.8 to 1.4 M . Stellar ages were estimated 
by using the values of effective temperature and luminosity to place the stars in a grid of 



3 ht tp : / / ams .astro . univie . ac . at / ~ vald / 
4 A(Fe)= log [N(Fe)/N(H)]+12 

5 Available at |http: / / verdi.as.utexas.e du/moog.html 
6 http://kurucz. harvard.edu/ 
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isochrones from Girardi et al. (2000) for [Fe/H] = 0.00 and 0.20, with ages from 0.063 to 
7 Gyr. In addition, stellar masses were calculated using the relation between luminosity, 
surface gravity, and effective temperature: 

M L g T e g 
log XT = log ~j— + log 4 log- , (1) 

M & L Q g Q Jeff,© 



with T e fr i0 = 5777 K and log g = 4.44. The masses and ages estimated for the targets stars 
are presented in Table [31 
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Table 1. Observing log. 
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V 
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N 
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Date 
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Table 2. Selected Fe lines and Measured Equivalent Widths. 
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7189.146 


Fc I 


3.07 


-2.771 


40.4 


35.3 


53.7 


37.5 


45.5 


38.5 


62.8 


6084.111 


Fc II 


3.20 


-3.780 


22.0 


32.2 


18.5 


31.7 


32.8 


21.6 


31.0 


6113.322 


Fc II 


3.22 


-4.110 


12.5 


21.0 


9.6 


19.1 


19.8 


12.0 


19.0 


7222.394 


Fell 


3.89 


-3.276 


19.7 






29.6 


29.1 


19.2 


27.1 


7224.487 


Fell 


3.89 


-3.226 


21.3 


32.4 


16.3 


31.8 


30.5 


20.0 
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3.2. Spectrum Synthesis 



Lithium isotopic ratios are measured via spectrum synthesis analysis. The Li I doublet 
at ~ 6708 A is an asymmetric blend of two 7 Li lines which consist of 7 hyperfine compo- 
nents, spanning a separation of approximately 0.15 A between stronger blue and weaker red 
components. The 6 Li lines present the exact same configuration, except with 3 hyperfine 
levels, with a total separation of approximately 0.15 A between stronger blue and weaker red 
components. As the 6 Li lines are much weaker, with the stronger one situated almost at the 
same wavelength as the weaker 7 Li component, the former isotope appears as a perturbation, 
providing an additional small asymmetry to the blend. In addition, the lithium region suffers 
from the contribution of several blends from metal and molecular CN lines. 



3.2.1. Broadening Parameters 

The stellar projected rotational velocity, v sin i, and macroturbulent velocity, V m , are 
broadening parameters which need to be defined for modelling spectral lines via spectrum 
synthesis. A good strategy to estimate rotation and macroturbulent velocities is to analyze 
lines which are isolated, unblended, and with similar strengths to the Li line and then apply 
the same broadening parameters to the Li I synthesis. A few Fe I lines which fall in the same 
echelle order as the Li I feature were investigated and the Fe I at 6703.567 A was identified 
as a clean line which was used to estimate v sin i and V m for the studied stars. 

A grid of synthetic spectra was computed for combinations of v sin i and V m varying 
between and 10 km s -1 (with steps of 0.1 km s _1 ). Also, we let the iron abundance vary 
within 0.05 dex of the abundance value which was obtained for that line. Best fits between 
synthetic and observed line profiles were obtained from a x 2 minimization as follows: 

y 2 = 1 y-(°»-^) 2 (2) 

v ' i=i 

where 0« and Si are, respectively, the observed and synthetic normalized fluxes at a wave- 
length point i across the line profile; a is rms error of the continuum, given by (S/N) -1 ; 
d = n — p is the number of degrees of freedom in the fit, where n is the number of points in 
the observed spectra used in the fit and p is the number of free parameters in the calculation 
of the synthetic spectra. In this case, p = 5: v sin i and V m ; r (continuum level), w (wave- 
length) and iron abundance A(Fe). Small adjustments in the continuum level (r < 0.4%), to 
account for possible errors in the normalization process were allowed. In addition, shifts in 
the central wavelenghts of the Fe I lines were needed in order to properly match the observed 
lines. 
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Fig. 2. — Spectroscopic determination of the effective temperature and microturbulent ve- 
locity for HD 82943 obtained from zero slopes in the runs of Fe I abundances with excitation 
potential of the transitions (upper panel) and reduced equivalent widths (lower panel). Fe I 
(filled squares) and Fe II (open squares) abundances are consistent and the slopes are zero 
for A(Fe) = 7.79. 
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Since the resolution of the bHROS spectra is R ~ 143,000, the velocity resolution is 
just a bit larger than 2 km s -1 . Thus, at low values of v sin i and typical macroturblent 
velocities in stars of the type studied here, the spectra will not be sensitive to changes in 
low projected rotational velocities found in the program stars. In order to estimate a lower 
limit to v sin i that can be detected significantly with the data here, synthetic spectra were 
computed for the Fe I 6703.567 A line with a range of values of v sin i, from 1.0 to 3.0 km 
s _1 , a single macroturbulent velocity of 3.0 km s" 1 and additionally smoothed to the bHROS 
resolution, using a model with solar parameters. The synthetic spectra were then sampled 
at the same scale as bHROS, after which noise was added such that S/N = 1,000. These 
"degraded" synthetic spectra were then subjected to the same analysis as the real spectra, 
with the result that for v sin % values of 3.0 and 2.5 km s _1 the proper rotational velocities 
were recovered to within 0.3 km s _1 (close to what the estimated uncertainties are for the 
target stars). For values of 2.0 and 1.5 km s" 1 , the analysis yielded a lower and constant 
value of 1.3 km s -1 , indicating that the bHROS spectra are not necessarily sensitive to values 
of v sin i < 1.5-2.0 km s -1 for typical macroturbulent velocities. 

We note, however, that even though the bHROS spectra are not sensitive to projected 
rotational velocities of less than about 2.0 km s _1 , it is still found in the real data that small 
values of v sin % provide somewhat better fits to the Fe I line profiles than a combination 
of larger v sin i and lower macrotubulence. Although these lower values of v sin % provide 
better fits, their detection should not be considered significant, with realistic lower limits of 
v sin i = 1.5 km/s. 

Figure [3] shows the fit obtained for the star HD 74156. The top panel plot both the 
observed and best-fit synthetic Fe I line profile as an illustration of the quality of the fits that 
are obtained in the analysis. The bottom panel illustrates the ^-minimization techniques 
used for estimating both projected rotational velocity (left) and macroturbulence (right). 

The uncertainties in the overall line broadening caused by macroturbulent and projected 
rotational velocities can be estimated by varying v sin i (while keeping V m fixed) in an interval 
of 1.0 km s _1 around its best value, with steps of 0.05 km s -1 . For each value of v sin i, Ax? 
— Xr " Xrmin was calculated and the velocity change which produced a la change in Ax? 
was taken as the uncertainty. The uncertainties in V m were obtained in a similar manner. 
The adopted broadening parameters and uncertainties are presented in Table HJ 
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Table 3. Adopted Stellar Parameters and Derived Metallicities. 



Star 


^eff 


log g 




M 


Age 


A(Fe) 




(K) 


(cm s -2 ) 


(km s -1 ) 


(Mq) 


(Gyr) 




HD 17051 


6197 


4.49 


1.24 


1.2 


1.0 


7.73 


HD 36435 


5503 


4.56 


1.46 


1.0 


1.0 


7.52 


HD 74156 


6100 


4.36 


1.38 


1.4 


2.5 


7.64 


HD 82943 


6055 


4.56 


1.33 


1.2 


1.5 


7.79 


HD 147513 


5904 


4.63 


1.48 


1.0 


2.0 


7.52 


HD 217107 


5690 


4.44 


1.25 


1.1 


5.5 


7.87 



Table 4. Elemental abundances and isotopic ratios. 



Star 


v sin i 


v m 




A(Li) 


6 Li/ 7 


Li 




(km s _1 ) 


(km s 










Sun 


1.70 


2.80 


0.96 


0.00 


HD 17051 


4.90 ± 0.25 


5.25 ± 


0.30 


2.48 ± 0.01 


0.03 ± 


0.04 


HD 36435 


4.25 ± 0.20 


2.95 ± 


0.30 


1.60 ± 0.03 


0.06 ± 


0.08 


HD 74156 


3.10 ± 0.25 


4.85 ± 


0.25 


2.59 ± 0.01 


0.00 ± 


0.03 


HD 82943 


< 1.5 


3.60 ± 


0.10 


2.49 ± 0.01 


0.00 ± 


0.02 


HD 147513 


< 1.5 


3.05 ± 


0.10 


2.03 ± 0.01 


0.02 ± 


0.03 


HD 217107 


< 1.5 


2.85 ± 


0.10 


< 0.36 


0.00 ± 


0.04 
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Fig. 3. — Top Panel: Observed (filled squares) and synthetic (solid line) spectra for the 
selected Fe I line which was used to define the broadening parameters in target star HD 
74156. The Xr-minimization used for estimating v sin i (filled circles; left) and V m (filled 
squares; right) is shown in the bottom panel. 
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3.2.2. Convection and line asymmetries 



The additional small radial velocity shifts noted above represent convective s hifts which 

are re lated to the effects of granulation in the stellar atmospheres (see discussion in lAUende Prieto et al 
20021 ). Figure H] shows the trend of radial velocity shifts applied as a function of the mea- 
sured equivalent widths for target star HD 82943 and the Sun. The trend of increasingly 
positive radial veloc it ies a s the equivalen t widt hs inc rease has been noted previously by 
Allende Prieto et all J2002h . iReddv et all J2002h . and iMandell et all J2004T ). All three of 
these studies measured the Sun and fit straight lines to the trends with slopes and inter- 
cepts respectively of: 2. 9 m s~ l (mA)" 1 an d 80 m s _1 ( Allende Prieto et al.ll2002); 2.8 m s -1 



mA)- 1 and 120 m s" 1 (IReddv et all 120021 ): 4.2 m s" 1 (mA)" 1 and 30 m s" 1 (IMandell et al. 



2004J). Since different Fe lines were used by the different groups, as well as different fitting 
techniques, the differences between the trends are not surprising and do not lead to signifi- 
cant differences in convective shifts. The results reported here for the Sun are a slope of 2.0 
m s _1 (mA) -1 and an intercept of —80 m s _1 . 

Solar-type stars are known to dis play slight red asymmetrie s in spectral lines due to the 
convective motions of granules (e.g. lAllende Prieto et al.ll2002l ); the flux across the line is 
dominated by hot, rising granules, while the cooler, falling (i.e., red-shifted) inter-granule 
regions produce the small red asymmetries. The line-bisector for the Fe I 6703.567 A line 
is shown in Figure 0, where it is plotted as flux level versus velocity (in m s _1 ) instead of 
wavelength. A detailed discussion of stellar line bisectors can be found in Gray (2005). The 
filled squares are the observed points in HD 82943, while the open squares were derived 
from a synthesis of this region using a standard ID model, which necessarily produces a 
symmetric line profile. The synthetic profile shows a vertical bisector, indicating an iso- 
lated symmetric line (demonstrating that the Fe I 6703.567 A line is good for determining 
broadening parameters). The observed line-bisector in HD 82943 deviates to the red due to 
convective granules; the magnitude of the convective asymmetry is manifested as a ~ 200 
m s _1 excursion in the line bisector. The shape and amplitude of this bisector is typical for 
stars of this type and represents a rather small perturbation of the line profile. 



3.3. The Line List 

The determination of the 6 Li/ 7 Li isotopic ratio rests upon modeling the shape of the 
Li I line-profile, which depends not only on the isotopic ratios themselves, but also on the 
stellar broadening mechanisms discussed in the previous section. In addition, the shape of 
the Li I profile is affected by several weak absorption features that fall within the wavelength 
region of the various lithium transitions. There are several weak lines of CN and other metal 
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Fig. 4. — Convective shifts measured for Fe I lines in the solar spectrum (top panel) and HD 
82943 (bottom panel). The solid lines in each panel represent the least-squares fits to the 
points. 
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Fig. 5. — Synthetic (open squares) and observed (filled squares) line bisectors measured 
for Fe I at 6703.567 A. The asymmetry to the red in the observed profile is indicative of 
photospheric velocity fields. Errors in the bisectors are largest in the central portion of the 
line. 
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absorption lines that will blend with 6 Li, thus a detailed compilation and careful assessment 
of line lists is a critical ingredient for the analysis. 

As mentioned in the introduction, the search for 6 Li in planet hosting stars has been 
the subject of a few recent studies and line lists for the Li I region are available in the 
literature. As a starting point i n this analysis , we ad opted the line list from the most recent 
and complete study to date by iMandell et al.l (12004 ). The list of lines and atomic data (<?/- 
val ues and wavelengths) carefully compiled by the authors were adjusted in that study to fit 
the iKurucz et al.l (11984] ) solar spectrum. In addition, a laboratory carbon arc spectrum was 
used to adjust the CN lines (in both wavelengths and (//-values) in comparison to theoretical 
line lists. As for the 3 possibilities for the unidentified feature at 6708.025 A listed by Mandell 
et al. (2004; Si I, Ti I and Ti II) we adopted the Si I line. Recall that Mandell et al. (2004) 
found no significant differences in the 6 Li/ 7 Li ratios derived using the 3 different lines. The 
first step here was to return all atomic data in their line list to those yalues appearing in th e 
original references. The CN molecular data were kept the same as in IMandell et al.l (120041 ) . 



This init i al line list was checked via comparisons between the solar flux spectrum from 
Kurucz et al.l (119841 ) and synthetic spectra which are based on the models and synthesis code 
adopted in this analysis. The solar spectrum was modelled using an ATLAS9 ODFNEW 
model with T e g = 5777 K, log g = 4.44, and £ = 1.24 km s _1 . The broadening parameters 
adopted were obtained from fits to the unblended Fe I line at 6703.567 A: v sin i = 1.70 km 
s \ V m = 2.80 km s -1 . The limb darkening coefficient, e, was taken from I Van Hammd (119931 ). 
however the choice of e has little effect on the Li I line profile and no measurable impact 
on the estimation of 6 Li isotopic fractions. The comparison of observed solar spectrum and 
synthesis in the small wavelength interval around the Li I feature was improved by a small 
continuum adjustment of 0.25%. In addition, small wavelength corrections to the nearby Fe I 
line (—0.003 A) and to Li I (+0.002 A) line led to a better fit; these wavelength corrections 
are those expected due to convective shifts in the solar photosphere (see Section I3.2.2p . A 
small change in the Fe I (//-value of +0.015 dex improved the fit to the Fe I profile, while 
the accurate laboratory (//-values for the Li I components were maintained. 

Within the scope of the above changes, the fit to the solar spectrum exhibited some 
small mismatch, especially near 6708.275 A, where there app eared additiona l abso rption. 
The absorp tion in this region w as attributed to a Mg I line by IMandell et al.l (120041 ) and a 



V I line by iReddy et al.l ( 2002T). No referen ces to either of these lines could be found in a 



literature search, whereas King et al.l (119971 ) attributed the absorption at this wavelength to 
a Ca I line (present in VALD-2). There were other small differences between the line lists: 
in the narrow region from 6708.31 to 6708.54 A the Mandell et al. list contained 4 CN lines 
while King et al. had 2 (at 6708.375 and 6708.635 A). We included the Ca I line plus the two 
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adjacent CN lines; some of the 5/- values were changed somewhat and wavelengths allowed to 
shift by 0.001 A in order to improve the fit to the solar spectrum. The final line list adopted 
in this study is presented in Table EJ The best-fit to the solar spectrum achieved with this 
line list was excellent and this is shown in Figure O The best-fit was o btai ned for A (Li) 



= 0.9 6, which is very close to what was found by both lKing et al.l (1997) andlReddy et al. 
( 120021 ). and within 0.1 dex of the recommended value in lAsplund et al.l (120051 ). 
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Table 5. The Line List 



A 


Identification 


LEP 


log 


;gf 


Original 


(A) 




(eV) 


(d 


ex) 


Reference 


6706.548 


CN g 2 (93) (11,5) 


3. 


.130 


-1 


.359 


M04 


6706.567 


CN g 2 (80) (8,3) 


2. 


490 


-1 


.650 


M04 


6706.657 


CN R l2 (22) (7,3) 


0. 


.870 


-3 


.001 


M04 


6706.733 


CN Q 1 (22) (7,3) 


0. 


.870 


-1 


.807 


M04 


6706.844 


CN i2i(34) (12,7) 


1 


.960 


-2. 


.775 


M04 


6706.863 


CN P 2 (83) (7,2) 


2 


,070 


-1 


.882 


M04 


6706.880 


Fe II 


5. 


.956 


-4. 


.504 


V 


6706.980 


Si I 


5 


,954 


-2 


.797 


V 


6707.205 


CN g 2 (47) (11,6) 


1 


,970 


-1 


.222 


M04 


6707.282 


CN g 2 (60) (10,5) 


2. 


.040 


-1 


.333 


M04 


6707.371 


CN Qx(85) (12,6) 


3 


.050 


-0 


.522 


M04 


6707.431 


Fe I 


4. 


.608 


-2. 


.288 


R02 


6707.457 


CN P 12 (13) (7,3) 


0. 


.790 


-3 


.055 


M04 


6707.470 


CN Qi(28) (12,7) 


1 


.880 


-1 


.451 


M04 


6707.473 


Sm II 





,933 


-1 


.477 


V 


6707.518 


VI 


2. 


,743 


-1 


.995 


V 


6707.545 


CN g 2 (44) (6,2) 





,960 


-1 


.548 


M04 


6707.595 


CN g 2 (29) (12,7) 


1 


,890 


-1 


.851 


M04 


6707.596 


Cr I 


4 


,208 


-2 


.767 


V 


6707.645 


CN P 21 (44) (6,2) 


0. 


,960 


-2. 


.460 


M04 


6707.740 


Ce II 


0. 


.500 


-3 


.810 


R02 


6707.752 


Scl 


4 


.049 


-2. 


.672 


V 


6707.756 


7 Li 


0. 


.000 


-0 


.428 


H99 


6707.768 


7 Li 





.000 


-0 


.206 


H99 


6707.771 


Cal 


5. 


,796 


-4. 


.015 


R02 


6707.807 


CN i2i(64) (5,1) 


1 


.210 


-1 


.853 


M04 


6707.848 


CN i?!(61) (19,12) 


3 


,600 


-2. 


.417 


M04 


6707.899 


CN P 2 (39) (20,13) 


3 


,360 


-3 


.110 


M04 


6707.907 


7 Li 





.000 


-1 


.509 


H99 
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4. LITHIUM ABUNDANCES AND ISOTOPIC RATIOS 

As discussed previously, several blending lines due to metals and CN affect the shape of 
the Li I feature in metal-rich stars. It is possible to study which of these potentially offending 
lines have the largest effects on the 6 Li absorption via test spectrum syntheses which can 
isolate each of the nearby lines. 

Synthetic spectra were generated to focus on lithium and each of the possible blending 
species in turn: CN, Si I, Ca I, and V I. In the case of CN, the main impact on 6 Li comes 
from 6707.807 A which falls slightly blueward (~ 0.10 A) of most of the 6 Li absorption and 
does not affect the derived abundances significantly. The Ca I line at 6708.275 A falls too 
far to the red to affect 6 Li, while V I at 6708.094 A is too weak to even be detectable. The 
biggest effect on 6 Li absorption comes from the blending line due to Si I at 6708.023 A, where 
its absorption falls on the red side of possible 6 Li absorption. Even here, absorption from 
6 Li that is about 1% deep would compete with Si I (at solar metallicity) and be detectable. 
These tests indicate that 6 Li absorption which is ~ 1% deep would be marginally detected, 
while larger amounts of 6 Li that would produce absorption of only a few percent would be 
detectable. 

The important goal is to obtain an overall good fit across the entire Li region from 
6707.3 A to 6708.4 A and this includes fitting the nearby feature just blueward of Li I which 
is mainly Fe I with some contribution from CN. Also, a close region nearly free of spectral 
lines (6706.40 - 6706.55 A) is useful for providing an estimate of the local continuum level, 
while the wavelengths of Fe I (6707.431 A) and Li I were adjusted slightly to take into 
account convective shifts (see Figure H|). The abundances of Ca and V were allowed to 
vary within ±0.2 dex in order to improve the overall fits to the observed spectra, as well as 
adjusting C and N to fit the blue CN feature. The final values, or limits, to 6 Li were set by 
a x^-minimization of A(Li), 6 Li/ 7 Li and A(Si) simultaneously. 

Lithium abundances and isotopic ratios are given in Table H] for the six target stars. 
The uncertainties estimated for the Li abundances and isotopic ratios were derived in a 
manner similar to that used for the rotational and macroturbulent velocities. The total 
lithium abundance was varied (but keeping 6 Li/ 7 Li fixed) in steps of 0.001 dex, with the xl 
computed for each value. A change in Xr of 1 was taken to define the ±lcr- value of A(Li). 
The uncertainty in the isotopic ratio was determined in an analogous way, with 6 Li/ 7 Li 
values varied in steps of 0.005 and the total Li abundance being fixed. The errors in A(Li) 
and 6 Li/ 7 Li are also shown in Table HI 

In order to further test the sensitivity of our method to the isotopic ratios, synthetic 
spectra were computed for 6 values of the 6 Li/ 7 Li ratio (between 0.00 to 0.05 with a step of 



Table 5 — Continued 



A 


Identification 






Original 


( A ^ 
( A ) 




(eV) 


(dexj 


Kelerence 


6707.908 


7 Li 


0.000 


-0.807 


H99 


6707.919 


7 Li 


0.000 


-0.807 


H99 


6707 920 


6 Li 


0.000 


-0.479 


H99 


6707.920 


7 Li 


0.000 


-0.807 


H99 


6707 923 


6 Li 


0.000 


-0.178 


H99 


6707.930 


CN CW35) (12,7) 


1.980 


-1.651 


M04 


6707.964 


Til 


1.879 


-6.903 


V 


6707.980 


CN i2 2 i(72) (10,5) 


2.390 


-2.027 


M04 


6708.023 


Si I 


6.000 


-2.910 


103 


6708.026 


CN P 2 (35) (12,7) 


1.980 


-2.031 


M04 


6708.073 


6 Li 


0.000 


-0.303 


H99 


6708.094 


VI 


1.218 


-3.113 


V 


6708.147 


CN P 2 (42) (11,6) 


1.870 


-1.434 


M04 


6708.275 


Cal 


2.710 


-3.377 


K97 


6708.375 


CN 


1.979 


-1.097 


K97 


6708.499 


CN 


1.868 


-1.423 


K97 


6708.577 


Fel 


5.446 


-2.728 


V 


6708.635 


CN P 2 (42) (11,6) 


1.870 


-1.584 


M04 
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Fig. 6. — Synthetic solar spectrum (solid line) computed with the line list in Table [51 The 
main lines cont ributing to th e synt hesis are identified. The observed solar spectrum (filled 
squares) is from lKurucz et al.l (11984J ) . The bottom panel shows the differences between model 
and observations; the overall agreement is excellent. 
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0.01) adopting the stellar parameters derived for HD 82943. These model spectra were then 
sampled at the same scale as bHROS, after which noise was added such that S/N = 750 
(typical S/N for our observations). These "degraded" synthetic spectra were then subjected 
to the same analysis as the real spectra and we were able to recover isotopic ratios within 
±0.01, demonstrating an analysis sensitivity of 0.01 for 6 Li/ 7 Li. 



As an additional test, we used the line list in llsraelian et al.l ( 120031 ) in order to analyze 
th e target star HP 82 943. The line list adopted in that study is identical to the one used 
by iReddy et al.l (120021 ). except for having a Si I line at 6708.025 A instead of a Ti I line. 
All the input parameters (model atmospheres, broadening parameters and convective shifts 
for the Li I and the nearby Fe I line) and the analysis method employed were the same 
as described above; the only difference being the removal of the Ca abundance as a free 
parameter because the feature at ~ 6708.275 A is attributed to V I in Israelian et al.'s list. 
The best fit between model and observations was achieved for A(Li) = 2.49 and 6 Li/ 7 Li = 
0.01, which are in excellent agreement with the obtained results using the line list in Table [5j 
T he derived isotopic rat io for this star, however, is significantly lower than the value derived 



in 



Israelian et al.l (120031 ) : although the Li abundances in the two studies agree well. 



5. DISCUSSION 

5.1. Three-Dimensional Model Atmospheres and Their Impact on the 

One-Dimensional 6 Li Results 



As discussed in Section 13.2.21 convective granules result in small red asymmetries in 
the spectral lines of solar-type stars, such as those analyzed here for 6 Li. Since the 6 Li 
components of the neutral lithium ground-state doublet fall to the red of the 7 Li transitions 
and because 6 Li is expected to be much less abundant than 7 Li, the presence of 6 Li will 
result in a slight red asymmetry to the combined Li I feature. Since both convection and 6 Li 
can result in a red asymmetric profile, a discussion of the expected convective asymmetries 
to the Li I feature must be included in order to interpret either detections or upper limits to 
6 Li fractions. 

A more realistic theoretical treatment of the convective motions that cause red-asymmetric 
line profiles has been included in m odel stellar atmospheres v ia the incorporat i on of radia- 
tive hydrodynamical convection, e.g. IStein fe Nordlundl (119981 ) or lAsplund et al.l (120001 ) . The 
resulting model atmospheric structures are three-dimensional, time-dependent and contain 
self-consistent convective flows; such models, often referred to as "3D models", naturally 
produce absorption lines with small red asymmetries. Recent analyses of the Li I line pro- 
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Fig. 7. — Observed and synthetic Li I profiles for HD 82943. The syntheses were calculated 
for 6 Li/ 7 Li isotopic ratios of 0, 5, and 10 %. The Chi-squared minimization of the 6 Li/ 7 Li 
isotopic ratio is shown in the bottom panel. The best fit is obtained for a synthetic spectrum 
with zero contribution from 6 Li. 
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files in stars that span the T e ff-range of the program stars here (although mostly for much 
more m etal-poor halo dwarfs and subgiants) have included a discussion of effects due to 3D 
models (Usplund et alihood : Cayrel et al. 2007,2008). 



Asplund et al.l (120061 ) determined 6 Li abundances for a set of 24 metal-poor halo dwarfs 



and subgiants based on very high quality spectra acquired with the ESO VLT/UVES. The 
authors found detections of 6 Li at the 2a significance level in 9 stars. They also showed 
that the isotopic ratio 6 Li/ 7 Li is effectively immune to changes in the stellar parameters 
(at the levels of ±100 K, ±0.2 dex, ±0.2 dex, and ±0.5 km/s in T eff , log g, [Fe/H], and 
microturbulent velocity, respectively). Another important result is that the usage of 3D 
model atmospheres in LTE (when compared to ID models also in LTE) does not change 
significantly the derived 6 Li fractions from synthesis of the Li I 6707.8 A line, suggesting 
that ID models are essentially as good as 3D models when employing LTE analyses. Finally, 
the authors observe that ID 6 Li/ 7 Li isotopic ratios are, in general, similar or lower than the 
ones derived from 3D LTE calculations. The average 3D-1D difference is ±0.008 and the 
largest value is ±0.033. In only 3 the 3D values smaller (—0.010, —0.015, and 

—0.017) than the corresponding ones for the ID case. These results suggest that the upper 
limits for the isotopic ratios set here would probably increase if the analysis had been done 
with 3D. However, the maximum increase would be compatible with the typical errors in 
the isotopic ratios. 



More recently, ICayrel et al.l (120071 ) showed that convection-induced line asymmetries 
could mimic the presence of 6 Li, at th e few percent level, for the halo turn-off metal-poor 
star HD 74000. Note that ISmith et all Jl998h found 6 Li/ 7 Li = 0.00±0.02 for this star via a 



ID analysis, similar to the analysis performed in this study. The ICayrel et al.l (120071 ) result 
was b ased on a 3D hydr odynamical simulation for a single set of atmospheric parameters, 
while ICayrel et al.l (120081 ) presented theoretical line asymmetries for a ra nge of atmospheric 
parameters that cover the values of the stars from Asplund et al. ( 2006h. Note tha t a 3D - 
NLTE treatment, the most adequate for the Li I analysis according to lAsplund et al.l (120061 ). 
is used by Cayrel (2007, 2008). The authors comp ute an asymmetry of about 2%, which is 
exactl y the mean value of 6 Li/ 7 Li for the stars from I Asplund et al. (l2006h . Also, ICayrel et al.l 



(120081 ) investigate the dependence of these asymmetries on the stellar parameters T eff , log g, 
and [Fe/H]. Based on a small grid of 8 hydrodynamical simulations, they predict an increase 
in the Li I line asymmetry with metallicity for [Fe/H] > —2.0. Also, they find that hotter 
stars with lower gravities show larger asymmetries than cooler stars with higher gravities. 

Giv en the tight li mits on 6 Li fractions derived here via ID measurements, combined 
with the ICayrel et al.l (120081 ) predictions that there may be expected line asymmetries of the 
order of or larger than ~2% for more metal-rich stars, there is no evidence of detectable 6 Li 
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in the 6 stars studied here. The average value of 6 Li/ 7 Li for the 6 stars in Table 5, inversily 
weighted by their respective uncertainties, is 6 Li/ 7 Li = 0.012, which is consistent with zero 
6 Li within the errors. If the non-planet hosting star, HD 36435, is rejected from the average 
due to its low S/N spectrum (Tabled]), the weighted average for the 5 stars with planets is 
even lower: 6 Li/ 7 Li = 0.008. With limits close to zero, these results also can provide tests 
on the 3D models and NLTE analyses. 



5.2. Limits on Accreted Mass 



An approximate upper limit to the number of 6 Li atoms in the atmospheres of the 
target stars can be estimated by combining limits to the number of 6 Li atoms with the 
convective zone masses. To carry out this exercise w e use t he convective zone masses for 
main sequence stars presented in iPinsonneault et al.l (120011 ) in their Figure 1. This is a 
simple estimate which assigns a single convective zone mass to the target stars (see table E]) 
based on their effective temperatures, but will provide a rough upper limit on the number 
of 6 Li atoms present in the outer layers of stellar atmospheres. The 6 Li abundance is set by 
the 6 Li/ 7 Li values in Table HI in those stars where this value is 0.00, the limit is set to the la 
uncertainty. Table [6] presents the resulting limits to the number of 6 Li atoms. In addition, 
the 6 Li abundances are translated to upper limits of accreted mass in Jovian masses. This 
is an admittedly naive estimate, but does provide a framework in which to limit the amount 
of accretion that could have occurred on the surfaces of these 5 planet-hosting stars. 

Accretion limits of a fraction of a Jovi an mass are typical values derived by both 
Mandell et al.l (120041 ) and iReddy et al.l (120021 ) . The tightest accretion limit set here is for 
the coolest star, HD 217107, due to the increasing strength of the Li I line with decreasing 
T e fj; this star has a very low total lithium abundance with a corresponding low limit of 6 Li. 
The tighter limit on HD 217107 is perhaps also the most interesting from this sample, as 



Table 6. Mass Accretion Limits 



Star 


log M C z 


N 


M acc 




(Ms) 


(*Li) 


(M Jup ) 


HD 17051 


-2.44 


<4.6xl0 43 


<0.33 


HD 74156 


-2.23 


<1.0xl0 44 


<0.71 


HD 82943 


-2.15 


<6.2xl0 43 


<0.44 


HD 147513 


-1.90 


<3.8xl0 43 


<0.27 


HD 217107 


-1.66 


<2.9xl0 42 


<0.02 
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this star is the only one with a very closely orbiting massive planet, with a sin i = 0.073 AU 
and M p i anct sin i = 1.33 M Jup . 

While the non-detections of 6 Li, at levels 1-2% of total lithium, are secure limits, the 
interpretation of this lack o f 6 Li in light of accretion rests on details of stellar evolution. 



Montalban fc Rebold (120021 ) have presented calculations of the expected evolution of the 



(i 



Li abundance as a function of time after ingestion for various stellar model masses and 
metallicities. Within the standard model framework of stellar evolution, where there is no 
exchange of material between the convective envelope and the radiative interior, accreted 
6 Li would survive for gigayears in mai n-sequence stars with M > 1.0-1.1 M and having 



metallicities of [Fe/H] = 0.0 and +0.3 (IMontalban fc Rebold 120021 ) . These types of stellar 



masses encompass most of the target stars studied here, thus the limits would suggest that 
any accretion of planetary material was typically less than a few to several tenths of a Jovian 
mass. 

When stellar evolution with non-standard transport processes, such as microscopic dif- 
fusion or turbulent mixing, is used to investigate the fate of 6 Li on the surface of solar-typ e 



stars, the interpretation becomes more complex, as discussed by IMontalban fc Rebold (120021 ). 
Both diffusion and turbulent mixing act to move material between the convective surface 
layer and the deeper radiative interior, with the result that accreted 6 Li will be removed from 
the photosphere and destroyed. The efficiency of this destruction is a function of stellar mass, 
metallicity, and the time at which 6 Li is accreted. For warmer effective temperatures, T c g > 
6100 K, the survival of detectable quantities of 6 Li can exceed a Gyr, thus if accretion takes 
place after such stars settle onto the main sequence, there would be measurable amounts of 
6 Li in some of these stars. Large, sensitive surveys for 6 Li remain useful ways to probe and 
constrain accretion. 



6. CONCLUSIONS 



The main conclusion drawn from the analysis presented here is that no detections of 6 Li 
are found in 5 planet-hosting stars which span the T e g range of 5700-6100 K, masses from 
~1.1— 1.4 M , and metallicities from [Fe/H] ~ +0.1 — 1-0.4. Since the hotter, more massive 
of the solar-type stars have the least massive convection zones, the stars studied here would 
be some of the best candidates for detecting signatures of accretion. The combination of 
high spectral-resolution with high-S/N makes this s earch one of the most sensitive and the 
results here can be combined with previous studies (IReddy et al.l 120021 ; iMandell et al.l 120041 ) 
to provide strong limits on accretion of 6 Li. 
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